NASA/TP-2017-219597 


Neutron Angular Scatter Effects in 
3DHZETRN: Quasi-Elastic 


John W. Wilson 
Old Dominion University, Norfolk, Virginia 


Charles M. Werneth, Tony C. Slaba 
Langley Research Center, Hampton, Virginia 


Francis F. Badavi 
Old Dominion University, Norfolk, Virginia 


Brandon D. Reddell 
Johnson Space Center, Houston, TX 


Amir A. Bahadori 
Kansas State University, Manhattan, KS 


—EpERhmn aEEaE—E—EE——EE— EEE —EEEE>e>>——————SSaSa-__a_a_a_a_a_a_a_a_a_a_a_a_a_a_aa_aa_—_E__—__==_ 
March 2017 


NASA STI Program .. . in Profile 


Since its founding, NASA has been dedicated to the 
advancement of aeronautics and space science. The 
NASA scientific and technical information (STI) 
program plays a key part in helping NASA maintain 
this important role. 


The NASA STI program operates under the auspices 
of the Agency Chief Information Officer. It collects, 
organizes, provides for archiving, and disseminates 


NASA’s STI. The NASA STI program provides access 


to the NTRS Registered and its public interface, the 
NASA Technical Reports Server, thus providing one 
of the largest collections of aeronautical and space 


science ST] in the world. Results are published in both 


non-NASA channels and by NASA in the NASA STI 
Report Series, which includes the following report 


types: 


e TECHNICAL PUBLICATION. Reports of 


completed research or a major significant phase of 


research that present the results of NASA 


Programs and include extensive data or theoretical 


analysis. Includes compilations of significant 
scientific and technical data and information 
deemed to be of continuing reference value. 
NASA counter-part of peer-reviewed formal 
professional papers but has less stringent 
limitations on manuscript length and extent of 
graphic presentations. 


e TECHNICAL MEMORANDUM. 
Scientific and technical findings that are 
preliminary or of specialized interest, 
e.g., quick release reports, working 
papers, and bibliographies that contain minimal 
annotation. Does not contain extensive analysis. 


e CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


e CONFERENCE PUBLICATION. 
Collected papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or co-sponsored by NASA. 


e SPECIAL PUBLICATION. Scientific, 
technical, or historical information from NASA 
programs, projects, and missions, often 
concemed with subjects having substantial 
public interest. 


e TECHNICAL TRANSLATION. 
English-language translations of foreign 
scientific and technical material pertinent to 
NASA’s mission. 


Specialized services also include organizing 
and publishing research results, distributing 
specialized research announcements and feeds, 
providing information desk and personal search 
support, and enabling data exchange services. 


For more information about the NASA STI program, 
see the following: 


e Access the NASA STI program home page at 
http://www.sti.nasa.gov 


e E-mail your question to help@sti.nasa.gov 


e Phone the NASA STI Information Desk at 
757-864-9658 


e §=6Write to: 
NASA STI Information Desk 
Mail Stop 148 
NASA Langley Research Center 
Hampton, VA 23681-2199 


NASA/TP-2017-219597 


Neutron Angular Scatter Effects in 
3DHZETRN: Quasi-Elastic 


John W. Wilson 
Old Dominion University, Norfolk, Virginia 


Charles M. Werneth, Tony C. Slaba 
Langley Research Center, Hampton, Virginia 


Francis F. Badavi 
Old Dominion University, Norfolk, Virginia 


Brandon D. Reddell 
Johnson Space Center, Houston, TX 


Amir A. Bahadori 
Kansas State University, Manhattan, KS 


National Aeronautics and 
Space Administration 


Langley Research Center 
Hampton, Virginia 23681-2199 


SSS eee 
March 2017 


The use of trademarks or names of manufacturers in this report is for accurate reporting and does not 
constitute an official endorsement, either expressed or implied, of such products or manufacturers by the 
National Aeronautics and Space Administration. 


Available from: 


NASA STI Program / Mail Stop 148 
NASA Langley Research Center 
Hampton, VA 23681-2199 
Fax: 757-864-6500 


Contents 


ADSULACT ss cases stbesiseh eves eh tinted ti dle teioe Sie ioh eusens ie aa ao shee nooo ete Sad Nettie aaa Lats eee 1 
Tintin Ct Ons see tecs de ssssss steed couch So Sapa ctes ev sestebeancbands sass ds sep bus Lovkss Sedsgnstbei gh aaeds shay ates sveuas do yagnsdbves oat so bpecbas (oshes betbapsoeedssaselesiees 1 
Determumistic:Code: Development es /osscsscecsacsctesscesseesieesehcesserdesssedseg soeecse at vevsepdesaua tein select eg eadheeteeseb conden aes setoeteene nies 1 
Simple Homogenous Spherical Geometry... cee eee esse cseecseeeeeeeeeeeeeeeceeccesecesecesecsaecsaecsaecsaecaeseaeseeeeeeeesesseeaeesaeeeaees 8 
A Dynamical Theory of Quasi-Elastic Scattering 0... eee ceecesecsecseeceecseeeeeeeeeeseesseeeseesecaecsaecsaecsaecaeseaeseaeeeeeeaeeees 10 
New: Quas1-Elastic- Model. ans HZET RN iiss ccseeccdtc bs estes cucehe cheat cck cha chu w Suete deh etob ba feubee ahs Cates Dee dhe dubes Cod an oeatee bavevacsdeetooes 17 
COmCLUSIONS 5 sees ce-ssssehsa sess d coach Sep sbis Lessa bets ze sted sacs de Sbsesbus JovbabsSasncteedsbaah da sepeates Sei gabde bageadbedsoosh bo bea odes Gysasbetbanedeeisess age 18 
ACkniowled Sements: i, cdcsaess edesd deus eecoescasseyastesb cusses ev vevsvtcusgeateevitdebeseessdeatveesepdhsstcrcgi eavech sanueedveu use duper dues ydtestiuees Ane 18 
REfETONCeS:, estes MAE ee aces sR SR aed a aes ao aE Neha a oh a Nas ade Sed hated NN as ee Me aah ee 18 


ill 


Figures 


1. Isotropic and forward fraction of neutrons produced in nuclear reactions by nucleons incident on aluminum 


according to the current HZETRN database... eee eeeeesseceseeesssecesneecsaeceneecnaeceneecsaeceeeecsaeceeneecnaeceeneesaeeeeees 4 
2. Isotropic scattering and forward scattering fractions of neutrons in aluminum according to the current elastic 
database‘and’ ENDF database... :s..ssscscs.s: sscueessseeseeus os casssbestesassopsosesnedeasasdupchpasoteduat siupsnee done Sebatey chp setevsegaeensassaaas 4 
3. “Geometry for 3D transport procedure... :....:..ssccscesce.cesessseessesnte vase sssenseistey cb eusevetvensoupenensosdonvasbepcbpasssevesesvepsoeeseds 6 
4. Spherical geometry and external source orientation for benchmark ComparisoNs..........0...:esseeeeeeeeeeeeeeeeeeeees 9 
5. Neutron fluence spectra at 35 g/cm? (left) and 40 g/cm? (right) in sphere exposed to Webber SPE. Results 
from 3DHZETRN (N = 18) using the for/iso and ge/mp models are shown along with MC evaluations. ....... 9 
6. Neutron production cross sections for 100 MeV proton collisions in the ge/mp and for/iso formalism. ...... 10 
7. Transmitted quasi-elastic spectra of a multiply scattered nucleon at four depths in nuclear matter as 
evaluated by equations (37) and (44) at four incident energies. 0... eee eeeceeneeceeeceeececeseceeaeeceaeceeaeecsaeeeeneers 14 
8. Impact of current Serber model of quasi-elastic contribution to neutron fluence spectra in an aluminum 
SPHOLe sores eda Sssieast ose tis ee esaeete god sddeesush va cu eag ie clue, Sun atiacs edema uceets cette daa Oks clgae Sia davies cupesas dee eh cease pede gue Mace aye 17 
9. Comparison of equation (3) quasi-elastic model with the new quasi-elastic model... eee eee eee estes 17 


iv 


Tables 


1. Root mean square nuclear charge radius, a- (fm) [Hofstadter and Collard 1967]. 0.0.0... ceesceeesceeneeeeteeeeneees 15 


Abstract 


The current 3DHZETRN code has a detailed three dimensional (3D) treatment of neutron 
transport based on a forward/isotropic assumption and has been compared to Monte Carlo (MC) 
simulation codes in various geometries. In most cases, it has been found that 3DHZETRN agrees 
with the MC codes to the extent they agree with each other. However, a recent study of neutron 
leakage from finite geometries revealed that further improvements to the 3DHZETRN formalism 
are needed. In the present report, angular scattering corrections to the neutron fluence are 
provided in an attempt to improve fluence estimates from a uniform sphere. It is found that further 
developments in the nuclear production models are required to fully evaluate the impact of 
transport model updates. A model for the quasi-elastic neutron production spectra is therefore 
developed and implemented into 3DHZETRN. 


Introduction 


The approach taken throughout the development of the HZETRN radiation transport code has been to 
develop a progression of transport solutions from the simple to increasingly complex, allowing early implementation 
of high-performance computational procedures based on marching algorithms [Wilson 1977; Wilson et al. 1986a, 
1991, 1994; Slaba et al. 2010a] and establishing a converging sequence of approximations with well-defined 
accuracy criteria. In support of this paradigm, a preliminary nuclear database was established in the early 
development efforts (circa 1973), so that numerical solutions of the Boltzmann equation could be investigated and 
developed. This early work is recounted in Wilson et al. [1988a,b; 1991]. The preliminary nuclear database used for 
numerical transport solution studies proved to be reasonably adequate for application to space system design in the 
sense that uncertainty was mainly associated with transport solution procedures rather than the nuclear database. 
Notably, the nucleon portion of the database as described by Wilson et al. [1991] has remained in the HZETRN code 
(and carried into 3DHZETRN) without any updates. Verification and validation processes of the nuclear and 
transport models were accomplished over the last several decades using realistic model solutions, convergence 
criteria, Monte Carlo (MC) simulation, laboratory experiments, and flight data [Wilson et al. 1991, 2005, 2007; 
Slaba et al. 2013; Matthiaé 2016]. Reviews of this development process, including enhanced three dimensional (3D) 
neutron transport, is given by Wilson et al. [2014a-c, 2015a-c]. 

Limitations of the current 3DHZETRN code and nuclear database were revealed in part by a study of 
neutron leakage in finite objects [Wilson et al. 2015b,c] wherein differences in leakage from cube and sphere 
geometries are obscured by the forward/isotropic (for/iso) approximation used in 3DHZETRN [Wilson et al. 
2014a,b]. Although 3DHZETRN and various MC codes are still plagued by uncertainty in nuclear models and data 
[Wilson 2014a-c, 2015a], the next step in code development appears as much limited by nuclear model uncertainty 
as by the transport methods. Therefore, specific problems identified in the neutron leakage study [Wilson et al. 
2015b,c] are targeted herein, and some nuclear model improvements are initiated. 

In the present report, the for/iso approximation in 3DHZETRN is replaced by a quasi-elastic/multiple- 
production (ge/mp) separation that allows angular dependence in the neutron production to be easily carried into 
existing transport procedures. Only minor modifications to the transport formalism and numerical methods are 
required to implement the ge/mp model. It is found that the quasi-elastic component of the neutron production cross 
section carried over from the preliminary nuclear model of Wilson et al. [1991] is inadequate for this application, 
and a more detailed model for quasi-elastic nucleon production is considered as a replacement. The new quasi- 
elastic model is described and implemented into 3DHZETRN transport procedures, and it is shown that modest 
improvements in the neutron fluence estimates for spherical geometry test-cases are achieved. 


Deterministic Code Development 


The relevant transport equations are the coupled linear Boltzmann equations derived on the basis of 
conservation principles [Wilson 1977; Wilson et al. 1991, 2005] for the differential flux (or fluence) density 


(x,Q,E) of aj type particle in the continuous slowing down approximation (CSDA) in which atomic processes are 


described by the stopping power S,(E) for each ion type j (vanishes for neutrons j = 7) as 


1a ae 
OV A, ag HP) + oF) 6, 6,(x,Q, B) =e. fon (E, E',0,2')6,(a,2',E')dQ'dE', (1) 


and solved subject to a boundary condition over the enclosure of the solution domain. In equation (1), o(£) is the 
total macroscopic cross section for a type j particle with kinetic energy E, and O it (E, £',Q,Q") is the double 


differential macroscopic production cross section for nuclear interactions in which a type k particle with kinetic 
energy E' and direction 92' produces a type j particle with kinetic energy F and direction 2. 

The double differential cross sections in equation (1) may be generally separated into reactive (r) and 
elastic (el) components according to 


O(B,E',0,2') = 09 (B, B',9,2') + cE, B',0,0'). (2) 


In equation (2), the elastic component is null except for k = j and describes a nuclear collision in which the projectile 
scatters from the target, but neither the projectile nor target break apart or become excited as a result of the 
interaction. The reactive component includes all other non-elastic nuclear processes including projectile and target 
fragmentation, de-excitation, and evaporation. The total cross section in equation (1), oj(E), may also be separated 
into elastic and reactive components as in equation (2). 

For transport code development efforts leading to 3DHZETRN, the reactive component of the differential 
cross sections for neutron production (j = n) was further separated as [Wilson et al. 2014a,b] 


nk,iso 


o) (E, B',9,0') = 0? ,,,(B, E',9,0") + ao, (E,E'). (3) 


The forward component of the cross section (for) is associated mainly with higher energy processes and projectile- 
like secondaries and is peaked near the projectile energy, E’, and direction, Q'. The isotropic component (iso) is 
mainly associated with lower energy particles produced, including evaporation products, but extends up to 
moderately high energies. Separation of the cross section into forward and isotropic components was carried out 
using the Ranft angular factor [Ranft 1980] and will be described later in this report. 

An alternative, but related, separation of the reactive cross section for neutron production is given as 
[Wilson et al. 1988a, 1991] 


o)(B B',.0,0') = 0) _(BE',9,0) +09 (BE',9,9). (4) 


Onk “qe Onk, ,mp 


In equation (4), the quasi-elastic term (ge) represents the spectrum of scattered particles having originated from the 
projectile and suffering at least one intra-nuclear collision with a nucleon within the target. Charge exchange 
processes are included in the quasi-elastic term. The multiple-production term (mp) accounts for those escaping 
nucleons and other light ions resulting from intra-nuclear collisions of the quasi-elastic scattered primary particles 
with nuclear matter and are associated with lower-energy particles produced including target fragments. 

Although similar in form, equations (3) and (4) bear important distinctions. Equation (3) provides a 
convenient separation of phase space into a high energy forward directed component and a lower energy isotropic 
component. The separation is achieved by identifying backward production represented in the Ranft angular factor 
as half of the isotropic spectrum. Subtraction and further manipulation then allows the forward component to be 
deduced [Wilson et al. 2014a,b]. Details of the physical processes included in the forward and isotropic components 
are not explicitly identified but are assumed to be reasonably characterized within the Ranft angular separation. 

Conversely, equation (4) is a more fundamentally physics-based separation of the cross section where 
specific physical processes and models may be assigned to the quasi-elastic and multiple-production components. 
One immediate advantage of equation (4) is that angular dependence in the lower energy multiple-production term 
can be maintained and carried into transport procedures. This work is focused on extending 3DHZETRN to utilize 
equation (4) instead of equation (3), thereby enabling a more detailed description of angular factors for secondary 
neutrons. Related efforts are then carried out to improve upon the quasi-elastic component of equation (4). 

In past versions of HZETRN and the current version of 3DHZETRN, the quasi-elastic component of 
equation (4) for nucleons is represented with the parametric form [Wilson et al. 1988b] 


N ae (E") 
C ve" +e 200-£/B') F (5) 


9n(0, E, Ap )o Oe E') 


of (E,B',0,0) = 97(0,E, Ap)oy” 


where Cy = 1 — In(2) / 20 is a normalization constant, NV, ,(£') is the quasi-elastic multiplicity (i.e. average 


jk,qe 
number of type j particles produced as a result of a quasi-elastic interaction) [Wilson et al. 1991], and of) (E') is 


the reactive component of the total cross section appearing on the left hand side of equation (1). The Ranft angular 
factor is given by [Ranft 1980] 


Npexp(-0? / Ap) ,0<0< 1/2 


° 6 
Npexp(—m / 4p) .7/2<0< 0 (6) 


(0, E, Ap) = 


where Nr is an energy dependent normalization factor defined such that J 9p (0, E, Ap )dQ = 1 and cos 6 = 2-0". 
The Ranft width factor is 


Ag = (0.12 + 0.000364) / E, (7) 


where Ar is the mass of the struck nucleus. 

The multiple-production term in equation (4) represents highly spectrally dispersed particles of lower 
energy including the de-excitation spectrum [Wilson and Costner 1977, Wilson et al. 1988b] with spectral 
components given by 


gl EEO. = 4184p) YY Ng BI en ‘ 
—n€ ne 


29, E ‘Aa Lien E ') 


(i) 
where NY re 


al. [1991]. 

As discussed previously in this section, early development efforts of 3DHZETRN [Wilson et al. 2014a-c, 
2015a], utilized the for/iso approximation to provide a convenient separation of phase space. Although not obvious 
in equations (3) - (8), the isotropic component of the cross section extends up to relatively high fragment energies 
and may include various physical processes with pronounced angular dependence. To illustrate this point, Fig. 1 
shows the forward and isotropic neutron production fractions for protons and neutrons incident on aluminum. The 
isotropic fraction is defined as 


E') are multiplicities from the Bertini model, and a; are shape parameters as described by Wilson et 


Es ee )(B, B',Q,0:') dE do 5 
me ee Jo (8, B'.9,2dEan 


where the integral in the numerator is only taken over the back hemisphere, 27B, and then doubled to represent the 
isotropic component of the cross section as defined by Wilson et al. [2014a-c, 2015a]. The forward component is 
simply one minus the isotropic component given in equation (9). Although the isotropic component dominates for 
projectile energies below ~20 MeV, as would be expected, it accounts for more than half of the cross section for 
projectile energies in the GeV region. Certainly, portions of the fragment spectrum for such high energy projectiles 
will include angular dependence that is not captured by the for/iso approximation. 


Solid lines: proton projectile 
Dashed lines: neutron projectile 
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Fig. 1. Isotropic and forward fraction of neutrons produced in nuclear reactions by nucleons incident on aluminum 
according to the current HZETRN database. 


A brief description of the elastic cross section is given here; although, no modifications to the transport or 
nuclear model for elastic interactions are considered in this work. The elastic scattering components are found using 
the differential angle distribution of Wilson et al. [1991] consisting of an S-wave amplitude and Chew's impulse 
approximation [Chew 1951] representing the higher order angular dispersion [Wilson et al. 1991] (similar to the 
approach of Wilson [1973] for the three nucleon problem). A linear energy factor is applied to the S-wave cross 
section to match the neutron KERMA as described in Wilson et al. [1991]. For further comparison, the fraction of S- 
wave scattering according to this model has been evaluated and compared with the Evaluated Nuclear Data File 
[ENDF 2016] derived S-wave fraction as shown in Fig. 2. It is clear that this simple scheme provides a reasonable 
first order approximation to the nuclear elastic scattering. 
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Dashed lines: ENDF 
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Fig. 2. Isotropic scattering and forward scattering fractions of neutrons in aluminum according to the current elastic 
database and ENDF database. 


Transport Model Development: Historical 


A prime limitation in developing solutions for equation (1) lies in evaluating the integral over dQ! at any 
arbitrary location within the media. The approach to a practical solution of equation (1) is to develop a progression 


of solution methods from the simple to increasingly complex allowing early implementation of high-performance 
computational procedures and establishing a converging sequence of approximations with well-defined accuracy 
criteria and means of verification and validation. The first step, leading to the lowest order solution, simplifies the 
evaluation over dQ! by introducing the straight-ahead approximation as guided by the nucleon transport studies of 
Alsmiller et al. [1965]. In these studies, MC methods were utilized, and the differential cross sections were 
approximated as 


04, (E, E',2,2') = o4(B, B)(Q-2'), (10) 


resulting in dose and dose equivalent per unit fluence to be within the statistical uncertainty of the MC result 
obtained using the fully angle dependent cross sections in slab geometry [Alsmiller et al. 1965, Wilson et al. 1991, 
2014a-c, 2015a]. Using asymptotic expansions about angular divergence parameters, it has been demonstrated that 
the error in using equation (10) for solving equation (1) is on the order of the square of the ratio of distance of 
divergence (few centimeters) to radius of curvature of the shield (few to several m) resulting in a small relative error 
in most human rated space systems [Wilson and Khandelwal 1974, Wilson et al. 1991]. The straight-ahead 
approximation of equation (10) formed the basis of the early versions of HZETRN and related numerical marching 
algorithms, for which the verification and validation processes are described elsewhere [Wilson et al. 2005, 2006, 
2014a-c, 2015a]. 

In the early Boltzmann transport studies of Wilson and Lamkin [1974, 1975] and Lamkin [1974], it was 
demonstrated that after neutrons are produced by a proton beam, the recoupling of neutrons back to the proton field 
is limited, but the neutron fields tend to build-up with increasing penetration depth. This mainly occurs since many 
protons produced by neutron-induced interactions are of lower energy and are quickly removed by atomic 
interactions (especially the isotropically produced protons). The probability of even a 100 MeV proton to undergo a 
nuclear reaction is only several percent [Wilson et al. 1991, 2005]. The low-energy neutrons on the other hand have 
no effective atomic interactions and propagate until a nuclear interaction occurs or they escape through the media 
boundary. Yet, the low-energy protons and other light ions produced in tissue through neutron collisions contribute 
significantly to biological injury and must be accounted for in transport methods [Foelsche et al. 1974]. This was the 
basis of the studies of Heinbockel and co-workers [2003] starting with Clowdsley et al. [2000, 2001] using the 
NUCFRGz?2 database [Wilson et al. 1987a,b; 1995; 1998]. 

Accordingly, improved methods for solving the neutron component of equation (1) were developed 
[Heinbockel et al. 2003; Slaba et al. 2010b]. These solutions approximated the neutron double differential cross 
sections as the sum of forward (f) and backward (b) components according to 


Onk (BE, E',9,') = Onk,f 


(E, B',9,2') + 0,4 (8, B',Q,2') . (11) 
Implementation of equation (11) resulted in the bi-directional code in common use for spacecraft design. A similar 
bi-directional approximation for the light ions was considered by Cucinotta et al. [1990] and Shavers et al. [1996] to 
study interface transition effects but was not integrated into the HZETRN software. 


Transport Model Development: Forward/Isotropic 
A sequence of approximations that span the formalisms from the simple straight-ahead and bi-directional 
methods to more complex propagation algorithms has been developed and is based on a simple for/iso assumption 


given previously by equation (3). The isotropic component is associated with lower-energy neutrons produced and 
defined by 


Rio EB') =2f o(B EA, dn, (12) 


where 27B represents the backward hemisphere. The forward component, 


ta r 1 ; is 
0? 4, (BE, B',0,0") = 0) (BE, E',A,0!) — reg. ) (EB, (13) 


nk,iso 


is associated mainly with direct quasi-elastic events [Wilson 1977, Wilson et al. 1988b], although other processes 
are included, and is highly peaked in the forward direction. 

The governing transport equation within the CSDA, ignoring multiple Coulomb scattering, is given by 
equation (1), and the total cross sections and double differential cross sections are the nuclear scattering and reaction 
values. The double differential cross sections are represented by equations (2) and (3). With the demonstrated 
success of the bi-directional approximation [Slaba et al. 2010b; Heinbockel et al. 2011; Wilson et al. 2014a-c, 


2015a] for neutron fields, the solution ; (a, 2, E) is divided into forward and isotropic components and the forward 


component is assumed to satisfy 


1 0 00 
QV — FSB) + (EB) 0) jor BB) = » J jor Bs EO, jor 9, B") AB", (14) 
: ; 


where the forward component of the cross section is given as 


0, (2, E',9,9') da + Jen 
Op for EB") = me : (15) 


* Jo efor B,B',0,2") + of(B, BY O,0)dQ > F= 0 


and the straight-ahead approximation is used to represent the angular spectra of equation (15). Note that QO only 
appears as a parameter in equation (14) and can be evaluated according to the boundary condition that is set as the 
incident fluence from direction Qo at the boundary as shown in Fig. 3. 


on 


74 


Fig. 3. Geometry for 3D transport procedure. 


The isotropic component of the neutron field is obtained by solving 


[RV + 0, (EN by igo (@ GE) = f° fon E,B' 2,20, inp, 2, B") dQ 4B 


= (16) 
oe is Onkiso Es E')y, fop (#2; Q%, E') dE! 
k 


using a bi-directional approximation along Q [Wilson et al. 2014a-c, 2015a-c]. Coupling between the forward 
solution of equation (14) and the isotropic neutron field is expressed through the perturbation given by the last term 


on the right hand side of equation (16). Although ¢ 


Sta (a, 2, E) is generated by the isotropic cross section term, the 
resulting fluence is not necessarily isotropic. A final step in the solution methodology is to compute the source of 
light ions produced from the lower energy isotropic neutrons. Once this source is computed, the isotropic component 
of the light ion flux is solved under the assumption that no further nuclear collisions occur, giving partial 3D 
treatment to low energy charged particles. 

The solution of these coupled equations (14) and (16) and some of their limitations are thoroughly 
discussed by Wilson et al. [2014a-c, 2015a-c]. The resulting transport algorithm allows direct evaluation of the 
straight-ahead and bi-directional approximation error at least in the context of the for/iso interaction approximation. 


6 


Although the for/iso formalism has yielded improved solutions over the straight-ahead and bi-directional models, 
limitations have also been discussed in previous sections and noted in recent work [Wilson et al. 2015b,c, Slaba et 
al. 2016] that will be addressed herein. The approach will be a more direct solution to equation (1) using an 
improved representation of neutron cross sections as in equations (4) - (7). 


Transport Model Development: Quasi-Elastic/Multiple-Production 
As a further refinement of detail for 3DHZETRN, equation (4) is now utilized to separate the reaction cross 


section for neutron production. In this more detailed approach, the fluence is separated into a high energy forward 
(hef) and angular (ang) component as 


?; (x, Q, BE)= Dj hep (@s 2 Ee) a Pj ang (® OQ, E) . (17) 


Equations (4) and (17) are substituted into equation (1), and the hef fluence is defined to satisfy 


1a 
SE pt a) Pete) =f 1. (E, Ey pep (0, E") dE", (18) 
J. 


where the cross section on the right hand side of equation (18) is written as 


ie oy, (E, E',2,0')d0 es 


ky (19) 
BR 0) (EEA) + MEE AM) ,j=7 


The angular spectra in equation (19) are represented within the straight-ahead approximation. The equation for the 
angular component of the neutron field is 


[Q+V +o, (E) Pn ang'®%E) = fo fon (Es E224, ang (@, O!, E') da" 
. < . (20) 
+f 0) ap (Es B', 2, %y) by neg (BE) AE" 
k 
Note that (@,Q,F) is generated by the angular dependent multiple-production interaction term, 


Pn, ang 
(EB, £',Q,.O) , and is not isotropic. 


The coupled equations (18) and (20) appear similar to equations (14) and (16) of the for/iso model with the 
exception that the cross section in the coupling source term is no longer isotropic. Despite this added complexity, 
equations (18) and (20) can be solved using the exact same numerical solution methodologies used in the for/iso 
model but with a modified source term. This will allow a direct evaluation of the errors associated with the previous 
forliso formalism. 

Equation (20) is solved by representing the double differential cross section of equation (20) within the bi- 
directional approximation, given by 


Onk “mp 


Onn (EB, B',9,.0' )= Onn, f (E, E')6(Q ~~ 2') + Onn pL, EO + 2') ? (21) 
Onn s B)E") =f. On (E, E'2,2)d0 = 0!) (EE) +0 p(B E) +o (BE), (22) 

and 
Onn» (Ey E') = fn? pom EB", 0,20 = 0) EB) + (BE), (23) 


where the forward (f) and backward (b) components of the multiple-production cross section are computed as 


OC mp) EE") = fo OO np (EB 8, QAO, (24) 
© J(BE)= (1) (B,E',A,2")dQ 25 
One taaiy ’ ) nB ort en , eee ) ‘ (25) 


In equations (22) - (25), 2nF and 27B represent integration over the forward and backward hemispheres relative to 
Q'. The subscripts (f) and (b) refer to the forward hemispheric and backward hemispheric components of the 
multiple-production neutron cross section, respectively. The forward and backward components of the elastic cross 
sections appearing in equations (22) and (23) are similarly evaluated by integrating the elastic values over the 
forward and backward hemispheres. 


Using equations (21) - (23) in equation (20) results in a coupled set of equations for Pn,ang (&s Q,£) and 
Dissaiig (az, —Q, E) . The equation for Drei (#, Q, E) is given by 
[RV +0, (E) by ang (®%E) = fn gE Ey ang (@ 2, E") dB 
+f Onn (Es Ey ang (@—O, B') aE" (26) 
yy ‘ cee (E, E , Q, %)Pr. neg (a, , E ') dE" 
k 
The equation for ¢,, ang (2, —Q, E) is 
[RV + 0, (BVI, ang(®—2,B) = fo ony Es Ey ang (®—O, B') dB 

+ Onno (Es En ang (; 2, B') dE" , (27) 


+)> is ae (£, E , —2, %) Pr nef (a, Q), E ') dE" 
k 


The second integral appearing on the right hand side of equations (26) and (27) is the coupling term between the Q 
and -Q multiple-production components. Equations (26) and (27) represent a coupled set of partial differential 
equations that are solved efficiently using a Neumann series approach as shown by Slaba et al. [2010b]. This series 


pee, 


is expected to converge quickly since 9, ai 


ie is required to conserve momentum. 

Similar to the for/iso formalism, a final step in the solution methodology is to compute the source of light 
ions produced from the neutrons of equations (26) and (27). Once this source is computed, the multiple-production 
component of the light ion flux is solved under the assumption that no further nuclear collisions occur, giving partial 


3D treatment to these charged particles. 
Simple Homogenous Spherical Geometry 


Simple geometry and source orientation are chosen to allow efficient MC simulation for comparison with 
with the above formalisms [Wilson et al. 2014a,b]. The geometry herein is taken as a simple aluminum sphere with 
a 40 g/cm? diameter. Target points are placed at depths along the z-axis. The external radiation environment is 
assumed to be anti-parallel with the z-axis, uniform in the x-y plane, positioned above the sphere, and directed down 
onto the top of the sphere as shown in Fig. 4. The February 1956 Solar Particle Event (SPE) spectrum [Webber 
1966] is assumed to be incident from above. Along with 3DHZETRN, the MC codes Geant4 [Agostinelli et al. 
2003], FLUKA [Fasso et al. 2005, Battistoni et al. 2007] and PHITS [Sato et al. 2013] were evaluated. Further 
details for the Monte Carlo codes and simulation setup can be found in Wilson et al. [2014a,b]. 
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Fig. 4. Spherical geometry and external source orientation for benchmark comparisons. 


The solution for N = 18 of both the for/iso and ge/mp formalisms are shown in Fig. 5 along with results 
from the MC simulations. Also shown is the hef component of the ge/mp solution to be compared with the for 
component of the for/iso solution. The most notable feature is the greatly increased total neutron fluence in the 10 
MeV to few hundred MeV domain in going from the for/iso approximation to the ge/mp approximation. The main 
difference appears to be in the hef and for components of the fluence. 

This difference arises for two reasons. First, the low energy angle distribution by Ranft becomes 
increasingly isotropic with decreasing fragment energy, and second, the quasi-elastic energy spectrum remains 
constant over most of its domain even at low energies as shown in Fig. 6 and equation (5). Therefore, the separation 
of the quasi-elastic component in the straight-ahead approximation is vastly different than the forward component of 
the for/iso approximation. The fault in the drop in the neutron fluence near 100 MeV in Fig. 5 results from the 
application of the for/iso approximation to the inadequate description of the nuclear induced spectrum, and an 
improved nuclear model is required. To begin correction of these factors we will develop an alternate and more 
reliable energy spectral distribution for the quasi-elastic cross section. 


10° 10° ; 
Geant4 FE Geant4 
FLUKA FLUKA 
10° ——— PHITS 10’ PHITS 
Reneterinetares for for” 
——_ for/iso A IVA Soriiso 
hef 6 eet PA AY hef 
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Fig. 5. Neutron fluence spectra at 35 g/cm? (left) and 40 g/cm? (right) in sphere exposed to Webber SPE. Results 
from 3DHZETRN (N = 18) using the for/iso and ge/mp models are shown along with MC evaluations. 
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Fig. 6. Neutron production cross sections for 100 MeV proton collisions in the ge/mp and for/iso formalisms. 


A Dynamical Theory of Quasi-Elastic Scattering 


In this section, a new model for estimating the quasi-elastic component of the neutron production cross 
section is described. Two solution methodologies are formulated to solve the model equations and are shown to give 
equivalent numerical results. A method for integrating the new quasi-elastic model into 3DHZETRN is also 
described. 

Historical high-energy nuclear reaction theory considered a rapid sequence of binary interactions within 
nuclear matter followed by a slower de-excitation process and allowed various reaction channels to be estimated 
[Serber 1947]. Initial implementation efforts of this model were based largely on semi-classical theory using mean 
free paths in nuclear matter and angular scattering models limited only by the Pauli Exclusion Principle (PEP) in the 
lowest energy and momentum transfer processes. The modern versions of this development are the intra-nuclear 
cascade (INC) codes utilizing computer generated random number sequences (initial attempts were two dimensional 
simulations using random number tables [Goldberger 1948]). Related to the INC model is a Boltzmann equation 
description, considered in this work, able to produce an analytic formalism equivalent to the INC stochastic results. 

The relevant transport equation is the linear Boltzmann equation derived on the basis of conservation 


principles for the flux density $(z,E) (nucleons/(fm?-MeV)) of nucleons moving through nuclear matter with kinetic 


energy E (MeV). Within the lowest order one dimensional approximation [Wilson et al. 1986b], the transport 
equation is 


0 ! ! ! 
2 + o(£)|¢(z,E£) = i o(E, E')¢(z, E")dE', (28) 


subject to the boundary condition ¢(0,£) = 6(E — E ) , with E,, being the lab energy (MeV) of the assumed mono- 
energetic beam with unit intensity incident on nuclear matter. Within the region of nuclear matter in which the target 
nucleus is represented, transport processes are modified by a liquid drop model with a potential well depth of Vo and 
a Fermi energy of & representing the top level of the Fermi sea below the binding potential, Vp, of the least bound 
nucleon. Note that the lowest energy transfers, A, given by A < &- &, are blocked by the PEP and that collisions on 
the range of ¢&- ¢ < A < Vo contribute to nuclear excitation. The solution domain of equation (28) is approximately 
over { & Ex-&}. For present purposes, the nuclear binding potential is ignored, and the integration limits are taken as 
{E, Ex-g } with g = 24.5 MeV. 

In equation (28), o(E) is the inverse mean free path in nuclear matter in units of 1/fm, and o(£,E’), in units 
of 1/(fm-MeV), represents the processes by which a nucleon with energy E' (MeV) collides with nuclear matter, and 


a nucleon emerges with energy E (MeV). The following parameterization is used in this work for o(E£), for E 
between 40 MeV and 500 MeV [Wilson et al. 1987b] 
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o(E) = E° / 16.6, (29) 


where the energy dependence is partially a result of PEP effects which reduce the cross section and increase the 
mean free path for a nucleon within the struck nucleus (especially at low energy), as has been determined 
empirically [Dymarz and Kohmura 1982]. For energies below 40 MeV and above 500 MeV, equation (29) is 
evaluated directly as an extrapolation. 

The cross section on the right hand side of equation (28), is further discussed in Wilson et al. [1986b, 1991] 
and approximated herein by 


BB!) = o(B)B(B) e+ o  B)B(E') 
o(E, E') = of E')B( ee + )B( rer (30) 
= 0, (£,E') +0,(£,E') 
where B(E’) is the slope parameter given by [Wilson et al. 1991] 
B(E') = 2me?10 3.5 4+ 30e 8! , (31) 


and mc? is the nucleon rest energy (938 MeV). The first term in equation (30), labeled as op(E,E’), accounts for 
elastic scattering of an incoming nucleon from a target nucleon and dominates for values of F near the incident 
energy E’. The second term, labeled as 0;(E,E"), provides the spectrum of recoil nucleons dislodged from nuclear 
matter and dominates for small values of FE. In the present formalism, details of nuclear binding effects (including 
Pauli blocking) and isospin are ignored for simplicity. Such corrections will be treated in a later work, once proof of 
principle is achieved in the present study. 

The operator on the left hand side of equation (28) can be inverted using an integrating factor and written as 
the Volterra integral equation 


z PE ‘ 
o(Z, E) = e 77 400, E) + J. f, 1 o(E)(z—2 )o(E, E'\(z"', E')dE' dz" (32) 


In order to establish the scattering properties of primary nucleons in nuclear media, solutions for equation 
(32) will be considered using only the oo(£,E’) component of the cross section in equation (30). It will then be 
shown in the following sections that such solutions may be used to compute improved quasi-elastic cross sections 
needed for 3DHZETRN, as discussed in the previous section. Equation (32) may be equivalently solved as a 
perturbation series [Wilson and Lamkin 1974, 1975; Wilson et al. 1986b] or with marching procedures [Wilson et 
al. 1988a]. The next two sub-sections describe both solution methodologies, and comparisons are provided as a 
verification of the numerical implementations. 

Note, the use of the term "multiple scattering" herein has a specific meaning and is more limited than that 
used elsewhere [Wilson 1974, 1975]. In this case, incident nucleons passing through the target may suffer one or 
more collisions with target nucleons; these multiple collisions are referred to as multiple scattering and are 
represented in solutions to equation (32) utilizing only the op(£,E’) component of the interaction cross section. 
Related to this terminology, the quasi-elastic cross section is defined in the next section in terms of the multiple 
scattering nucleon fluence transmitted through a specified finite region of nuclear matter. 


Perturbative Solution 
A perturbative solution methodology for equation (32) is developed in this section and is used to verify 
more efficient computational procedures based on a marching algorithm. The first term in the perturbative solution 


is obtained by neglecting the integrals on the right hand side of equation (32) and consists of the un-collided particle 
beam given as 
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dy (2, B) = e "60, E) 
e (4 §(E — E,) (33) 
00) (z,B) . 


Higher order terms in the series solution are found by successive iterations of the equation [Wilson and Lamkin 
1974, 1975; Wilson et al. 1986a] 


z pe, 1 ; 
E —o(E)(z—z') EF! (8) (41 B"\dE'dz' 
) = ie i, e€ y( ) Oe (2 ) )d dz', (34) 


for i> 1. The solution for primary nucleons suffering only one collision can be evaluated analytically and is given 
by 


(2, B) = is fore 0 (E, Ed?) (2', E\dE'dz' 
7 ~ -o(E)z _ 9-0 E,)2 (35) 
= (BB) o(E,) — o(E) 


The next term of the multiple scattering series represents those nucleons suffering two collisions and is given by 


=f oe (29 (B, B')d0(z', Bd" dz! 
et o,(E, E") Pl Ey)[eoe* — ero gol — gol) (36) 
— o( E') o(E,) — o(E) o(E") — of) 


The remaining integral in equation (36) is difficult to evaluate analytically but is easily handled with numerical 
quadrature. The remaining higher order terms are found by evaluating equation (34) numerically. 

Note that this is similar to the perturbation series solutions developed many years ago [Wilson and Lamkin 
1974, 1975, Lamkin 1974, Wilson et al. 1986b]. In the above, only the fate of incident nucleons is described, and the 
two nucleonic components {p,n} have been ignored. However, the formalism can be easily extended to include 
multiple constituents {p,n} and charge exchange cross sections (see Wilson et al. [1991]). The complete solution for 
the scattered primary nucleons that suffer at least one collision in nuclear matter is given as 


ins =D" z,E), (37) 


where the subscript ms is used to denote the multiple scattering perturbative solution. Note that the oy (z,£) term 
represents the incident fluence reaching depth z without suffering a nuclear elastic scattering event (neglecting the 
important diffractive components [Wilson 1974, 1975]) and is therefore not part of the multiple scattering solution 
in equation (37). 


Marching procedure 


Highly efficient numerical marching procedures have also been developed for equation (28). The solution 
proceeds by first noting that equation (32) can be equivalently written over a small step-size, h, given by 


bz +h, E) = ez, B) fe ie oy(E, Ez +t, E")dB' dt, (38) 
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where only the op component of the cross section is used, as in the previous section. Equation (38) propagates both 
the scattered and un-collided primary fluence over a step-size of h. However, for comparison to the multiple 
scattering solution of equation (37), the un-collided primary fluence needs to be removed from the marching 
procedure. This is achieved by noting that the fluence in equation (38) can be written as 


(2, E) = (2, E) + ¢,,(2,£), (39) 


where ¢0(z,E) represents the un-collided primary nucleons as in the previous section, and ¢;-(z,E) represents the 


nucleons suffering at least one nuclear collision. The solution for the scattered fluence, ¢,-(z,E), may be compared 


directly to the multiple scattering solution given by equation (37). The terminology (scattered vs. multiple 
scattering) and notation has been modified here to distinguish between the perturbative method and marching 
procedure solutions. 

If equation (39) is substituted into equation (38), one obtains 


Paez + A,B) = € PMG, (z, B) 

—o(E)h —o(E,)h 

+0) (E, E;, je (40) 
o(E,) — o(£) 


4 a ee: —o(E)(h—t) ! ' ' 
i ie e oy (E, E')d,,(z +t, E)dE dt, 


where equation (33) has been used to relate the un-collided primary nucleon field at an arbitrary depth to the 
boundary condition. An explicit marching algorithm is derived by noting that [Wilson et al. 1991] 


Pyol2 +t, B') = ONG, (2,.E') + Olt). (41) 
If equation (41) is substituted into the integrand of equation (40), one obtains 


se(2 + hy E) = eM". (2, E) 
~o(E)h —o(E,)h 
+0,(E, Eye 7 = = ee (42) 
o(E,,) — o(E) 
h E, ' 
wie, —o(E)(h—-t) ,—o(E')t ! ! ! 2 
J. : e e (ENG (E, Eb, (z, B)dE' dt + O(h?). 


SC 


Numerical requirements of using equation (41) are discussed in Wilson et al. [1991], where the marching procedure 
in equation (42) is shown to be unconditionally stable. As a result, the O(h’) notation is dropped hereafter. Over 
sufficiently small step-sizes, the integral from 0 to h can be approximated with a single-point expansion, which 
leaves the marching equation as 


-o(E)h _ ,—o(E,)h 
+04 (E, By )e-Ru)? a (43) 
o(E,,) — o(E) 
pe oak J “tg (E,E')@,.(z,B'"\dE' 
E 0 ’ sc\~) z 


Equation (43) is evaluated on a discrete energy grid with Nz points, allowing the integral to be approximately, but 
efficiently evaluated as 
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Puo(2 ar h, E) — eg (%, E) 
—o(E)h _ e EL )h 


+09(E, E,)e 7 © (44) 
ee o(E,) — o(B) 
—o(E)h Na 
theo) > (7), | Pscl% Ej) — Byo(% By 41)); 
j=i 
where 
= u Bisa ! ! 
(0), = gay Jn OE BME", (45) 
j+l j I 
®,.(,E) = [ by(2B "dB". (46) 


The solutions evaluated with the perturbative solution of equation (37) and the marching procedure of 
equation (44) are shown in Fig. 7 at four incident energies of Ez; = 100, 200, 300, 400 MeV and various depths in 
nuclear matter. It is clear that the most important component is near the incident energy at small penetration depths 
but diminishes at higher order scattering where further degradation of energy dominates at the larger depths. It can 
also be seen that the two solution methodologies give nearly identical results, although the marching procedure is 
more computationally efficient. 
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Fig. 7. Transmitted quasi-elastic spectra of a multiply scattered nucleon at four depths in nuclear matter as evaluated 
by equations (37) and (44) at four incident energies. 
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Connection to Ranft/Bertini Model in HZETRN 


Following Wilson et al. [1986b], the transmitted fluence spectra of scattered nucleons (equivalently 
obtained with either the perturbative method or marching procedure) can be utilized to estimate the differential 
energy cross section of quasi-elastically produced nucleons. Such estimates can then be integrated into 3DHZETRN 
to replace the currently used parametric form given by equation (5). Minor modifications to the Ranft/Bertini 
nucleon production model currently used in 3DHZETRN are needed to complete the integration process and are 
described in this section. 

Of particular importance is the evaluation of the absorption and the quasi-elastic cross sections for nucleons 
impacting target nuclei. Within the notation of Serber [1947], the absorption cross section (in units of fm) is related 
to attenuation of the scattered primary nucleon fluence incident on the target nucleus and given as 


o,(E,) = 2m Ap ele PL yay (47) 


where p is the target nuclear density function evaluated from the charge densities [Wilson 1975, Wilson and Costner 
1977, Wilson et al. 1991] normalized to unit density as accounted for by o(E£z) in units of fm‘. An equivalent 


definition of the absorption cross section, utilizing the un-collided primary nucleon fluence ¢0(z,E), understood as 


the probability of suffering no collisions in penetration to depth z, is given as 


Oqps(By) = 20 ie is J, Gol), Bidz bdb 


— On Ne “I — ¢ 7(Es x6) |bab 


(48) 


where ¢o was previously defined in equation (33), and z(b) is the penetration depth in a finite volume of nuclear 
matter for an impact parameter b. The absorption cross section defined in equation (47) or (48) includes both the 
reaction probability and average multiplicity and is consistent with the original definition for the absorption cross 
section given by the Serber model [1947]. 

Although equation (47) could be evaluated for any arbitrary density function, it will suffice herein to use a 
uniform sphere approximation allowing a simple evaluation as 


fe p(b,z)dz = VR? —b? = z(b), (49) 


where R = 1.29az, is the radius (fm) of the equivalent nuclear sphere [Wilson 1975], and the root mean square (RMS) 
radius, aa, is found from the RMS charge radii given in Table 1 using a, = ya — 0.64. Note that the function 


z(b) defined in equation (49) gives the chord length through a nuclear sphere for an impact distance b. 


Table 1. Root mean square nuclear charge radius, a, (fm) [Hofstadter and Collard 1967]. 
ac 
1 0.84 
2 217 
3 1.78 
4 1.63 
A 
2 


6<A<14 2.40 


16 0.58 + 0.82A"% 


Using equation (49), the absorption cross section is evaluated as 
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Say(E) = 2 f|1 - gre NEE bab 


we “IRI +4 20(E,)R|—1 > (50) 
TR? 4 


and the nucleon quasi-elastic differential cross section is found from 
a R 
6,(E,B,) = 20 | b,,[2(6),E]bab , (51) 


where ¢@,,, is computed for an incident beam with energy E, using either the perturbative method or marching 


procedure described in previous sections (marching procedure results are used hereafter for computational 
efficiency), and z(b) is given by equation (49). Note that equation (51) is the nucleon quasi-elastic cross section; the 
specific isotopic (i.e. neutron or proton) components are not yet represented. 

To make use of well-established nuclear reaction cross sections [Tripathi et al. 1997, 1998, 1999, 2002], the 
quasi-elastic differential cross section of equation (51) is re-normalized as 


o,.(E, E;) =o,( Ee) ; (52) 


where o;(Ez) is the nuclear reaction cross section from Tripathi et al. [1997, 1998, 1999, 2002] for a type k projectile 
with kinetic energy Ey. 

To preserve the isotopic distribution of the quasi-elastic cross section currently used in HZETRN, the 
fraction of type j particles (neutrons or protons) produced by collisions induced by type k projectiles (neutrons or 
protons) is evaluated as 


DY pe (E, ) 


Fieae\a) = $5, (E,) 
I 


(53) 


where Njxqe(Ex) is the quasi-elastic multiplicity currently used in HZETRN obtained from the balance of Bertini 
multiplicities as described by Wilson et al. [1991]. The quasi-elastic cross section for this new model is given as 


Op ge Er Er) = Fix ge Et )Oge(E, Ez) - (54) 


Note that in this new model, the ratio of quasi-elastically produced neutrons to protons remains the same as the 
current model in HZETRN [Wilson et al. 1991]; however, the overall magnitude and spectrum of the quasi-elastic 
component is determined by the new nuclear model. 

A final element of integrating this new nuclear model into 3DHZETRN requires the quasi-elastic 
multiplicity to be efficiently evaluated such that it can be re-balanced against the Bertini cascade multiplicities as 
described in Wilson et al. [1991]. The total quasi-elastic multiplicity is given by 


N ge (E) = Fie ge(#) (55) 


O abs (E) | 


where the term in brackets is the average number of primary nucleons transmitted through the target nuclear sphere. 

If this simple extension of the Ranft/Bertini model currently used in HZETRN is successful, then a more 
complete description could be at hand by extending the formalism to treat isotopic components (with Pauli blocking) 
and other particle production processes, as had been started many years ago [Wilson et al. 1986b]. The next section 
demonstrates the impact of this new formalism on the for/iso assumption and the ge/mp model. 
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New Quasi-Elastic Model in HZETRN 


The simple homogenous spherical geometry exposed to the Webber SPE considered in a prior section (see 
Fig. 4) is now re-evaluated with the new quasi-elastic cross section model. Using this new quasi-elastic cross section 
model in the for/iso and ge/mp transport formalisms of 3DHZETRN yields the results in Fig. 8. There are some 
encouraging improvements in the transport results with improved quasi-elastic cross sections in the 10 to few 
hundred MeV region. An expanded view of the prior quasi-elastic and the present quasi-elastic term is shown in Fig. 
9 where it is clear that this re-evaluation has improved the comparisons with the MC codes. Future work will 
evaluate the multiple-production cross section effects on the current solutions and will include Pauli blocking, 
affecting both the quasi-elastic and multiple-production processes. 
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Fig. 8. Impact of current Serber model of quasi-elastic contribution to neutron fluence spectra in an aluminum 
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Fig. 9. Comparison of equation (3) quasi-elastic model with the new quasi-elastic model. 
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Conclusions 


In this work the for/iso formalism in 3DHZETRN was replaced with a ge/mp formalism allowing angular 
dependence to be included in the multiple-production neutron source term. Using simplified spherical geometry 
exposed to the Webber SPE, it was found that the ge/mp transport formalism provides some improvement over the 
forliso assumption. However, further improvements to the quasi-elastic component of the neutron production cross 
section are needed. 

In order to improve the current quasi-elastic component, represented parametrically by equation (5), the 
Serber [1947] nuclear transport model was considered. Perturbative and numerical marching solutions were 
developed to solve the relevant equations representing transport of nucleons through a finite volume of nuclear 
matter. Details associated with nuclear binding were neglected. Methods for coupling the new quasi-elastic model 
into 3DHZETRN were also described. 

The spherical geometry benchmark case was then re-evaluated using the new quasi-elastic cross section 
model, and it was found that modest improvements were achieved. Further work will focus on refining the 
implementation of the Serber model to include effects associated with nuclear binding and Pauli blocking. 
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